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Abstract 

Yang-Mills theories undergo a deconfining phase transition at a crit- 
ical temperature. In lattice calculations the temporal Wilson loop and 
Z3 order parameter show above this temperature a behavior typical of 
deconfinement. A quantity of interest in the study of this transition is the 
gluon propagator and its evolution with temperature. This contribution 
describes the current status of an investigation of the finite temperature 
gluon propagator in Landau gauge. It analyzes the high temperature 
case. The resulting equations are compared to the corresponding ones 
of three-dimensional Yang-Mills theory. Under certain assumptions it is 
found that a kind of spatial "confinement" is still present, even at very 
high temperatures. 

1 Introduction 

The gluon propagator at zero and finite temperature is of many possible uses, 
although it is a gauge-dependent quantity. In Landau gauge the ghost propa- 
gator is of equal importance. Knowledge of these propagators can be used as 
input to phenomenological calculations, e.g. in analysis of heavy ion collisions 
and the properties of the produced state of matter. On the other hand these 
propagators can be used to compare analytical calculations with lattice sim- 
ulations to estimate the effect of approximations and to gain insight into the 
physical mechanisms underlying the lattice results and elucidate the dynamics 
of confinement. 

This contribution describes the approach and the status of the calculation 
of the finite-temperature gluon propagator in Landau gauge Yang-Mills (YM) 
theory. The gluon propagator in the energy regime of interest for heavy ion 
collisions is non-perturbative, and so must be the calculation. As infrared sin- 
gularities are anticipated, a continuum method is desirable. Therefore Dyson- 
Schwinger equations are employed as a non-perturbative approach. This will be 
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Figure 1: The truncated DSE used in the present work. The propagators with 
large full dots denote the fully dressed propagators while the ones without are 
the perturbative ones. A vertex denoted by a small full dot represents a pertur- 
bative vertex while the open circles indicate the fully dressed ones. In addition 
to the employed truncation all vertex functions are approximated by bare ones. 



described in section 2. The reason for choosing this method are the encouraging 
results at vanishing temperature which will be exemplified in section 3. Section 
4 describes the necessary extensions to access the equilibrium properties of glu- 
ons at finite temperature. Section 5 will deal with the high-temperature limit 
and discuss the comparison with a dimensionally reduced theory. Finally this 
contribution ends with a summary and an outlook. 



2 Gluon and ghost propagator equations 
2.1 Dyson-Schwinger equations 

Dyson-Schwinger equations (DSE) are a well-known tool for non-perturbative 
calculations in several areas of physics. The DSE for YM theories can be ob- 
tained by using functional derivatives H, ||. However, this approach produces 
an infinite hierarchy of equations: each equation of a given rt-point function 
depends on higher n-point functions. Since it is not possible to solve such an 
infinite coupled system, it is necessary to truncate it. There is no a-priori infor- 
mation about which truncation is sufficient to keep the relevant physics. The 
truncation scheme employed here is justified by comparison of the results at 
zero temperature with lattice results, which show that all qualitative features 
and many quantitative features survive. The truncation scheme is to keep only 
the equations for the 2-point functions and also only up to three-point vertex 
level as is depicted in figure |l|. 

In addition, it is assumed that multiplicative renormalization holds also 
in the non-perturbative regime and that therefore wavefunction renormaliza- 
tion and vertex renormalization can be cast into renormalization constants Zi. 
Renormalization is performed using a momentum subtraction scheme. The re- 
sulting equations after inserting the color structure and removing an overall 
factor of Sab are|] 

^^ik^D^^ik - q)G,Aq, k)DG{q) , (1) 
^AU quantities are defined using the conventions of ref. [^. 
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^^iq^Daik + q)DG{q)G,{k + g, q) 

Dai3{q)Dp„[k + q)Tp^p{~q,q + k, -k) . (2) 

To render the equations dimensionless, the dressing functions G and Z are 
defined via their relation to the propagators 

= D,Ak)^P,Ak)^ (3) 

for the ghost and the gluon, respectively, and the transverse projector is given 
by P^i, = k^kjj/k'^ . 

It remains to fix the fully dressed vertices, since they are unknown as long 
as the equations for the three point functions are discarded. It is possible to 
construct these vertices using the Slavnov- Taylor identities to constrain them as 
much as possible and this has been done At zero temperature, explicit calcu- 
lations show that vertices constructed in this way only slightly improve the result 
compared to perturbative vertices |^. The effect of other non-perturbatively 
dressed vertices has been studied in ref. |^. Since constructed vertices induce 
a significant complexity to the problem it is therefore more reasonable in this 
first approach to finite temperature to keep only the perturbative vertices. 



2.2 Gauge symmetry 

To obtain only scalar equations it is useful to contract the gluon equation with 
the projector P^,j. However, it turns out that such a contraction produces spu- 
rious quadratic divergencies. They stem from the fact that the truncation of 
the DSE violates gauge symmetry. To remove these divergencies and therefore 
the effects of gauge symmetry violation it is possible to use instead of the trans- 
verse projector other projectors which project onto states without gauge 
symmetry violations. This can be accomplished e.g. by using a generalized 
projector 

P',. - - (4) 

where C is a real parameter. Choosing C, = d = A removes the spurious quadratic 
divergencies . This procedure is not unique and introduces ambiguities in the 
value of numerical coefficients. This effect has been studied by varying the 
value of Q while removing manually the spurious divergencies. The results at 
zero temperature demonstrate that these effects are small |0| . 

A second item with respect to gauge symmetry in the non-perturbative 
regime are Gribov copies. It has been shown that for the solutions of DSE 
in Landau gauge it is sufficient to require positive semi-definiteness of the dress- 
ing functions G and Z for all momenta to stay within the first Gribov horizon 
1^. This is however not yet a full solution to the problem, since there are also 
Gribov copies within the first Gribov horizon. It is possible, in principle, to 
solve the problem by either introducing a second set of ghost fields to fully fix 
the gauge or by adding a new term to the DSE |^, ^. But, since lattice calcula- 
tions indicate that the infiuence of Gribov copies inside the first Gribov horizon 
is small ||lO|, these considerable complifications are neglected for now. 
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Figure 2: Comparison of lattice calculations |[T^ with results from DSE for 
Nc = 2. The left panel shows the dressing function of the gluon and the right 
panel the one for the ghost 

2.3 Kugo-Ojima confinement criterion 

When studying the finite-temperature gluon and ghost propagator one of the 
main goals to be achieved is the question of the fate of confinement. It is 
therefore necessary to find a criterion to test for the presence of confinement. 
Confinement means the absence of gluons from the physical spectrum. This is 
for example possible if they do not have a Kallen-Lehmann representation. Kugo 
and Ojima were able to construct a criterion to test for such a realization of 
confinement by assuming that the gluons carry BRST charge (l^. The criterion 
can be cast in a simple form in Landau gauge [|l2|: Does the eucHdean ghost 
propagator diverge stronger than a particle pole as the euclidean four momentum 
approaches zero, i.e. does G diverge for k ^ 0? To utilize this criterion and to 
simplify the calculations, the analysis will be performed in euclidean space. 

3 Results at zero temperature 

The calculations at zero temperature have been carried out in |^, ^, 0, Figure 
H shows a comparison of lattice calculations with results of the DSE. These are 
calculations for Nc = 2, but the dressing functions have the same form for 
Nc = 3, if 't Hooft scaHng is employed, since the DSE in this truncation order 
depend only on the combination g^Nc- 

Both propagators compare well with lattice results although the gluon 
propagator misses some strength around 1 GeV. It is likely that this is due to the 
neglect of the two- loop diagram in the equations for the two-point functions ||^. 
The results show confinement according to the Kugo-Ojima criterion with a 
divergence of G with an exponent of —k in fc^ with k = 0.5953. In addition, 
the gluon dressing function vanishes at zero momentum with an exponent of 2k, 
1^ 1^, ||. It is also possible to define a running coupling constant and it turns 
out that according to these calculations there is an infrared stable fixed point 
in YM theory. The good agreement with lattice calculations and the fact that 
the Kugo-Ojima criterion is satisfied motivates to use this method also at finite 
temperature. 
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4 Extension to finite temperature 



4.1 Formulation 

Since the primary interest are the equilibrium properties of gluons at finite tem- 
perature, it is convenient to use the Matsubara formalism. This induces that 
there are now two independent variables, the discrete zero component of the 
momentum and the absolute value of the three momentum. Additionally there 
are now two possible tensor structures for the gluon propagator ||l5| instead of 
only P^^. Both have to be four dimensional-transverse due to gauge symme- 
try, but one is three-dimensional longitudinal and the other three-dimensional 
transverse. They can be expressed as 



D — X k^ki, kijki, k^ko /i i 

^Tfii/ — O^u 72 ' "pO— 72 '""01^—72 f /jOCoi/ U + j 

3 3 3 3 

Plm. = Pf..~PTt.u (5) 



where k^ denotes the magnitude of the three momentum and the zero compo- 
nent of the momentum fco — 2TrnT is the bosonic Matsubara frequency. These 
projectors are orthogonal to each other and satisfy -Pr/i/i — 2 and Pl^^ = 1. 
There are therefore two independent dressing functions, Zl and Zt, and the 
gluon propagator is defined as 

P>f_,^[k) = Prtiiy—j^ h Phtiu ^2 ■ ("^ 

By this definition Zl and Zt have to become equal at zero temperature and 
equal to Z . 

It is possible to obtain two equations for Zl and Zt by contracting the gluon 
DSE with PT^lly and with PL^^ly respectively. After inserting the perturbative 
vertices, this results in the following set of equations for the three dressing 
functions: 



1 _ ~ 5 g'^TNc 

G(k) - ^^ + ^^72^ 

(AL(k, q)G{q)ZL(k -q) + AT{k, q)ZT(k - q)) (7) 



— Zt, — Zl'- 



J dddq^P{k,q)G{q)G{k + q) 



(NLik, q)ZL{q)ZLiq + k) + Ni{k, q)ZL{q)ZT{q + k) 
+N2{k, q)ZL{k + q)ZT{q) + NT{k, q)ZT{q)ZT{k + q)) (8) 

2^3-^1 ^^51 / dedq^R{k,q)G{q)G{k + q) 



ZT{k) 



{ML{,k,q)ZL{q)ZL{q + k) + Mi{k,q)ZL{q)ZT{q + k) 
+M2(fc, q)ZL{k + q)ZT{q) + MT(fc, q)ZT{q)ZT{k + q)) (9) 
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where T is the temperature and the sum runs over the bosonic Matsubara 
frequencies of the gluon and the ghost. The integral kernels Ai, P, R, Ni and 
Mi are quite lengthy and will not be quoted here. The trivial angular integration 
has been performed. 

Since no additional divergencies arise at finite temperature the renor- 
malization procedure can be kept by employing temperature- independent Zi. 

4.2 Gauge symmetry 

In the same sense as at zero temperature, spurious quadratic divergencies arise 
due to the truncation and projection of the DSE. It is possible to remove such 
artifacts in the same, and therefore also ambiguous, way as at zero temperature. 
Since there are now two independent functions for the gluon dressing functions, 
both projectors have to be modified and it is necessary to introduce two real 
parameters. This can be accomplished using the following generalized projectors 

Kg 

P^. = <:PT^u + {i-Ois^^-{^ + k)SM. (10) 

Kg 

It turns out that all spurious divergencies are removed for ^ = and C = 
3. These values are interesting, since they force the longitudinal projector to 
live only in the compactified dimension, while the transverse part becomes a 
structure which is reminiscent of the three-dimensional transverse projector. 

5 The high-temperature hmit 

The finite-temperature regime is quite complicated and it is therefore useful to 
start at some limiting cases and evolve their solutions to finite temperatures. 
One possibility is to start at zero temperature and to evolve the known solu- 
tion to non-zero temperatures. This approach is currently under investigation. 
The other approach is to start at high temperatures. As this method provided 
already some insights it will be described in the following. 

5.1 Definition of the high-temperature limit 

For the high-temperature limit, it is assumed that only the lowest Matsubara 
frequency contributes. It is therefore described by first letting the loop Mat- 
subara frequency go to zero and then the exterior frequency go to zero. The 
appearance of the explicit temperature dependence can be removed by either 
measuring the momentum in units of temperature or by defining a temperature 
dependent coupHng constant g^T — g|, which is dimensionful and its usefulness 
will be discussed in the next subsection. 

By performing this limit, several remarkable observations can be made. The 
first is, that many of the integral kernels vanish and that the complete system 
becomes independent of Their exact values are 
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o? sin^ 9 

Al = Q, At = , P = 0, 



Nl = 0, 7Vi = , N2 = , Nt = 



{ql - (C - l)ws -_C_ql cos2 9) sm9 



R= ' ' (11) 



^ - - 4.1 + 4(C - l)u,3 + 4C,i cos^ 9) sin. ^ 



sm f 



= ^fc|^((^3 + 2'?3)((C - 9)fc3' - ^4) + 8(C - 3)u;3(fc3^ + qi) 



I cos^ ( 



+ ((8(93 + (C + 7)fc^ + 4(5C - l)fc3'93) + 4(4C93' + ^KC + 3))u;3 + ^wl) 

where u± = A;| + g| ± 2^3930036' and W3 = fc3g3Cos6'. Secondly, the three- 
dimensional transverse projector now becomes the transverse projector of a 
theory with only three dimensions, while the longitudinal projector only acts in 
the compactified dimension. In addition, it is interesting to note that the value 
of C coincides with that of the three-dimensional Brown-Pennigton projector. 
But the most striking feature is, that if Z\ can be neglected compared to 
then the equation for Zl decouples. This leads to an interesting comparison to 
the three-dimensional theory. 

5.2 Comparison to three-dimensional Yang-Mills theory 

Compactified four-dimensional YM theory is not the same as three-dimensional 
YM theory, since the zeroth component of the four-dimensional gauge field be- 
comes an additional Higgs field However, lattice calculations indicate that 
this Higgs field is unimportant and produces only a small effect on the gluon 
propagator It is therefore interesting to compare to the three-dimensional 
YM theory only. The DSE can be calculated in the same way as in four dimen- 
sions and yield 

^ - " ' ( d9dqMh.qz)GM{qz)Z^d{h-q3) (12) 



G3<j(A:3) (27r 



2Z3 - Zi^^ J d9dq3P3ih, q3)G3diq3)G3d{h + Qs) 



+Zi0^ I d9dq3iN3ik3, (73)^3^(93)^3^(93 + h)) (13) 



with the integral kernels 



sm 

u2 2 



^3 = ,;^m + 2ql){kl{C~9)~4ql + 8w3{C-3){kl+ql) 



m^lt + (7 + Okt + AkiqiibQ - 1)) 

+Aws.{AC,ql + (C + 3)fc2) -f AC^wl) cos^ 9) (14) 
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using the three-dimensional Brown-Pennigton projector. It turns out, that the 
resulting coupled equations for the three-dimensional ghost and gluon are very 
similar. In fact, if Zj^ can be neglected, the three-dimensional coupling constant 
g\ is identified with g^T and the three-dimensional dressing function Z^d with 
Zt, then the above equations become identical to the ones for three-dimensional 
YM theory provided C = 3 is chosen. It is therefore tempting to identify the lon- 
gitudinal part of the four-dimensional gluon propagator with the contribution 
from the Higgs field. This is supported by the fact that the three-dimensional 
gluon propagator has also to be transverse due to gauge symmetry and the 
three-dimensional longitudinal part Zl therefore cannot be part of it. This in 
turn would justify the assumption that Zl is indeed small and can be neglected 
compared with Zt- However, only a detailed calculation of these quantities 
can prove this assumption. Additional support is found by the fact that in 
another approximation scheme, the Mandelstam approximation, the same situ- 
ation arises. Here is, however, no space to detail these other calculations which 
have been done by us. Moreover, recent lattice results support the finding of 
an infrared divergent ghost dressing function at finite temperature |Q which 
diverges weaker than at zero temperature as would be expected from the above, 
since in three dimensions k is smaller [||. 

If this turns out to be the case, it would have remarkable consequences, 
since the three-dimensional YM theory indeed does confine |^, and therefore 
a (spatial) kind of confinement would still be present even at very high tem- 
peratures in YM theory. This corresponds to lattice findings in that only the 
temporal Wilson loop shows deconfinement while the Wilson loops in the not 
compactified spatial dimensions do not ||l^ . This would underHne the fact that 
the magnetic sector of YM theory is non-trivial even at high temperatures. 



6 Summary and outlook 

This contribution describes our recent progress in a calculation of the finite- 
temperature gluon propagator in Landau-gauge YM theory. It is a first step 
generaHzing a self-consistent zero-temperature formulation within the Matsub- 
ara formalism. Possibilities to circumvent the problem of gauge symmetry viola- 
tions due to the truncation of the DSE have been discussed. A comparison of the 
high-temperature equations of the compactified theory with the corresponding 
equations for a lower-dimensional theory reveals interesting properties. If the 
assumption of a negligible infiuence of the three-dimensional longitudinal part 
of the gluon propagator at high temperature turns out to be correct, then this 
proves the non-triviality of the magnetic sector of the YM theory even at high 
temperatures. Note that this would imply the presence of some confinement 
effects at very high temperatures. The physical realization of such a spatial 
"confinement" has still to be understood. This may be possible as soon as the 
solutions of the DSE are obtained. The corresponding calculations are under 
way. 
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